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ABSTRACT 

We present the results of a near-infrared narrowband imaging survey of 
planetary nebulae. Objects were selected in a way that complements similar 
surveys done at visible and near-infrared wavelengths. No new detections of 
molecular hydrogen emission were made. The H 2 is frequently found to be 
extended, except in young, visibly compact objects. Our results are consistent 
with the already determined correlation of H 2 emission with planetary nebula 
morphological type. Filamentary and other kinds of structures are clearly 
resolved in many nebulae. 

Subject headings: ISM: molecules - planetary pebulae: general - stars: AGB and 
post-AGB - infrared: ISM: continuum - ISM: structure 

1. Introduction 

Ever since their discovery more than 200 years ago, planetary nebulae (PNe) have 
been the subject of intense study. These beautiful and enigmatic regions have revealed 
many of their secrets, most within the past 20 years. Recent reviews of PNe and PNe 
evolution include Kwok (1994), Peimbert (1992), and Pottasch (1992). With the advent 
of wide-field array detectors operating in the near-infrared (near-IR) spectral region, a 
new and potentially critical tool for the study of PNe is now available. While PNe emit 
throughout the electromagnetic spectrum, the near-IR has many advantages over other 
wavelengths. In particular, it is well-known that PNe contain a substantial quantity of dust 
that was probably formed while the star was on the asymptotic giant branch (AGB). The 
near-IR therefore provides lower attenuation to important dense regions than the visible. 
Atmospheric seeing at near-IR wavelengths is typically reduced compared to that at visible 
wavelengths. When coupled to high-spatial resolution cameras, this can help to resolve 
structure that is not normally seen at the shorter wavelengths. The molecular distribution 
in PNe is not well understood, and many transitions from vibrationally excited H 2 and 
other molecular species fall in the near-IR. While the H 2 transitions are excited by UV 
photons in the photodissociation front of some younger PNe (see Hora & Latter 1994; 
Dinerstein et al. 1988), for the most part H 2 in PNe is excited by shocks (Kastner et al. 
1994). Several studies have already been made that explore the near-IR properties of PNe, 
both as a group and of specific examples (e.g., Kelly & Latter 1995; Hora & Latter 1994; 
Kastner et al. 1994; Webster et al. 1988; Zuckerman & Gatley 1988). 
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It is for these reasons that we have undertaken the study to be discussed in this and 
following papers - a near-IR imaging and spectral study of a moderately sized sample of 
PNe. Our goals for this study were manyfold. Specifically, we have looked for systematic 
differences between the near-IR appearance of PNe to that in the visible (e.g., Schwarz et 
al. 1992; Balick 1987), as well as acquired new, high-spatial resolution data on extended 
H 2 emission. Because the near-IR transitions of H 2 are excited only in hot gas, or that 
exposed to intense UV fields, it probes the gas closest to the photodissociation region or 
that impacted by shocks. A focus of the study was on narrowband line imaging of H 2 and 
hydrogen Br7. Our object selection was made so that there is considerable intentional 
overlap with previous CCD and H 2 surveys. It has been suggested that there is a strong 
correlation of molecular detection (particularly H 2 ) with PNe morphological type, and 
possibly progenitor mass (Huggins & Healy 1989; Zuckerman & Gatley 1988). To more 
fully understand this correlation, attempts were made to make new detections of extended 
H 2 emission as well as image objects previously detected. A full interpretation of these data 
combined with the results of a near-IR spectral survey currently underway will be presented 
in additional papers. 

2. Observations 

The data presented in this paper include images acquired from Kitt Peak, Arizona and 
Mauna Kea, Hawaii between 1992 June and 1993 December. The Kitt Peak observations 
were made with the Steward Observatory (SO) 2.3m telescope using the SO 256x256 
HgCdTe NICM0S3 array camera and ~ 1% bandwidth tunable filters. The tilting 
narrowband filters were tuned to hydrogen Br7 (A = 2.1655 //m), H 2 u = 1 -> S(l) 
(A = 2.1213 /Wm), and adjacent continuum (see Latter et al. 1993 for details). Cold (77 K) 
reimaging optics provided an image scale of 0'.'6 pixel -1 . For some objects, broadband J, 
H } and K images were also obtained. Weather conditions during much of the observing 
period were highly irregular, making precise flux calibration difficult. Since determining 
the near-IR morphologies of PNe was a main goal, we have not flux calibrated the suspect 
data. In a typical 15 minute integration (on + off) in the narrowband filters, we found a 3<r 
sensitivity of 4.0 X 10 -16 ergs s -1 cm -2 pixel -1 in 0'.'6 pixels. Most integration times for 
images presented here were 6-10 minutes on source. For some objects such as J 900, NGC 
2346, and M 57, we integrated longer to bring out faint emission. The seeing at A = 2.2 /mi 
was in the range 0'.'7 - 1'.'5, with ~ l'.'l being typical. 

Observations on Mauna Kea were made using the University of Hawaii (UH) 2.2m 
telescope and UH 256x256 HgCdTe NICMOS3 array camera (Hodapp et al. 1992) at an 
image scale of 0'.'12 pixel -1 or 0'.'24 pixel -1 . Measurements were made in the J, H } and K 
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bandpasses, and fixed wavelength narrowband filters were used for observations of Br7, H 2 
v = 1 — > S(f), and A = 2.23 — 2.29 //m continuum. Because of the much broader bandpass 
of the UH continuum filter relative to the line filters, and possible H 2 line contamination, 
continuum subtraction was not attempted with the UH data. Observations were also made 
in the /-band using the UH 2048x2048 CCD camera. Seeing on Mauna Kea at A = 2.2 /mi 
was typically 0'.'5 - 0'.'7. Therefore, these data provide much higher spatial resolution than 
the Kitt Peak data (on average). For this reason, we have included data from both sites for 
some of the same objects. 

The images were constructed from many short exposures that were individually 
sky-subtracted and flat-fielded before being shifted and averaged to form final images. For 
extended sources, the telescope was nodded to nearby regions of blank sky to obtain sky 
frames. Compact sources could be observed efficiently by rastering them on the array. In 
this way, one object frame acted as the sky frame for the next. Flat-field frames were 
produced by averaging blank sky frames and then subtracting a dark frame of equal 
integration time. All of the data were reduced using standard IRAF routines. When 
accurate flux calibration was warranted, it was determined by comparison with standard 
stars from Elias et al. (1982). It is important to note that there can be a contribution from 
He I 3 3 P° — 4 3 S at A = 2.114 /mi in the H 2 narrowband images. In virtually all cases, such 
a contribution can be determined from the location of the emission. We will refer to this 
potential problem on a case-by-case basis. The data are summarized in Table 1, and filter 
parameters are listed in Table 2. The grayscale images are presented in Figures 1-15 
(Plates # - #)• Note that in the figure captions, UH data is specifically identified as such. 
Typically, four images are presented for each object selected from the following: /-band, 
J-band, //-band, A^-band, H 2 , Br7, narrowband continuum, continuum subtracted H 2 , 
and continuum subtracted Br7. 

3. Discussion 

We will break up our discussion into overall results and confirmed H 2 emission objects. 
Note that whenever possible, the morphological types listed in Table 1 and mentioned below 
are those of Balick (1987), or they follow the convention defined by him. A difficulty in 
morphological classification is that it can depend on the wavelength observed and flux level 
considered. In this paper, we will use the term "bipolar reflection nebula" to refer only to 
those post-AGB objects and proto-PNe (PPNe) that display a double-lobed appearance at 
visible wavelengths and have the characteristic that the lobes are seen primarily in reflected 
light from the central star (such as AFGL 2688). A "butterfly" or "bipolar" nebula is a 
PN that is bi-lobed and has an equatorial waist (e.g., M 2-9; see Schwarz et al. 1993). 
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These types are likely the same general structure viewed at different evolutionary states 
(such that a bipolar reflection nebula evolves to butterfly), but this relationship has not yet 
been proved. In general, we will use the term "polar lobes" to refer to those regions that 
are identifiable as the "wings" of butterfly-type nebulae or the lobes of the younger bipolar 
reflection nebulae. 

3.1. General Results 

The data in this paper represent a good selection of PNe morphological types. While 
many interesting features and structures have been found in a number of these objects, no 
clear systematic differences between visible light images and the near-IR are seen. This is 
meant in the sense that nebular classification in the visible is, in most cases, the same as 
that found in the near-IR. Differences seen are generally related to excitation of various 
species and the location of shocks. The Br7 emission follows Ha closely in objects where 
differential attenuation is not a factor, with apparent exceptions discussed below. The 
continuum in our ~ 1% wide filters can be significant, but typically the near-IR emission 
from PNe is dominated by Br7. For the most part, we have found that the near-IR 
continuum tends to mimic the Br7 distribution, and continuum subtraction does not offer 
significantly different morphological information. Given our sensitivity limits, extended 
halos of the type seen by Balick et al. (1992) could not be detected in Br7. 

Detailed structure is seen in all bands imaged. Bright knots or clumps can be found 
in most objects in both line and continuum emission. Filamentary structures have been 
resolved in a number of the PNe. Some PNe have relatively smooth halos surrounding 
complex inner structures. A transition from smooth, uniform mass loss to a much more 
dynamic process is indicated. The effects generated by the interaction of multiple velocity 
winds is clearly evident in some objects (e.g., NGC 7009). Others have structures that are 
difficult to unravel into a three-dimensional picture (e.g., NGC 6543, NGC 6826). In the 
following, we discuss several sources that were either not observed in H 2 or in which there 
is no clear evidence for H 2 emission. 

3.1.1. NGC 6302 

We have displayed images only of the central region of NGC 6302 (Fig. 5b), which is 
extended in the visible beyond the array limits (see Schwarz et al. 1992). However, data 
taken toward regions outside that shown indicate little extended emission, likely owing to 
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our sensitivity limits. Emission seen in the H 2 and the Br7 filters are very similar, while 
the continuum appears much less extended. This object has a very hot central star, and 
as such, it is unlikely that the emission observed in the 2.12 /mi filter is from H 2 , but is 
more likely He I 3 3 P° — 4 3 S contamination. The J-band image appears to have the same 
characteristics as the 2.26 fim continuum. However, a pronounced "arc" of emission is 
present on the south-east side, which does not appear so strongly in the other data. 

Lester & Dinerstein (1984) have presented maps of NGC 6302 in the J— and K— bands, 
together with longer wavelengths. They reported the detection of an infrared disk at the 
center of the nebula that might have been associated with the mechanism for collimating 
the outflow of material into the polar lobes. However, it is apparent from the higher 
resolution images presented here that the distribution of near-infrared emission is similar to 
the visible appearance, at least in the central region. In particular, the dark lane that is 
found in visible images is also seen in the J- and A^-bands. Therefore, we find no evidence 
for a disk or toroidal structure associated with the mechanism for shaping the nebula. 

3.1.2. NGC 6572 

NGC 6572 has been imaged previously in the visible (e.g., Jewitt et al. 1986; Schwarz 
et al. 1992). However, our near-IR images reveal much more detail (Fig. 6b). This object 
appears to be of the bipolar or early butterfly morphological type but is not compressed at 
the equator, so it is better to classify it an elliptical. There is an intensity enhancement in 
the equatorial region, however. Johnson & Jones (1991) found NGC 6572 to be polarized 
in the visible with P ~ 1.4%. This would suggest that at least part of the light observed in 
the near-IR is polarized as well. Indeed, there are indications of horn-like structures of the 
type seen in AFGL 2688 (see Fig. 14a). These features in AFGL 2688 have been modeled 
and found to be a result of single-photon scattering in an axisymmetric distribution of dust 
(Latter et al. 1993; see also Morris 1981). There is some evidence in the H 2 filter image for 
He I emission that appears as a bright spot to the northwest of the central star and is not 
seen in the other data. 

3.1.3. NGC 7026 

Emission at the location of the bright knots in our H 2 continuum subtracted images 
(Fig. 8b) is likely from He I. NGC 7026 appears similar in general morphology to NGC 
6572. Balick (1987) classified this object as a late elliptical. 
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3.1.4. BD+3CP3639 

This young, carbon-rich PN is similar in many respects to NGC 7027. Both are 
elliptical (in the infrared; Fig. 10) with strong near- and mid-IR excesses. Both have H 2 
emission in the halo regions, and the so-called unidentified infrared (UIR) emission features 
are prominent at 3.28 fim and in the mid-infrared (we did not image BD+30°3639 in H 2 ). 
BD+30°3639 has been imaged previously in the near-IR (Graham et al. 1993a; Hora et al. 
1993; Roche 1989; Smith et al. 1989). The A^-band image presented here is at higher spatial 
resolution and shows more details of the nebular structure. The nebula appears slightly 
rectangular, with the brightest portion of the lobes in the N and S regions. The central 
star is well separated from the ring at this resolution, showing clearly the evacuated central 
region. The upper right frame in Fig. 10b reveals a faint halo that extends out to about 
twice the diameter of the bright ring. The H 2 emission detected by Graham et al. (1993a) is 
visible at the E edge of the halo, as well as several other regions along the edge of the halo. 

3.1.5. IC 418 

IC 418 is a young, low-excitation PNe (Fig. lib) and has previously been imaged in 
the near-IR (Hora et al. 1993). Like BD+30°3639 and NGC 7027, it exhibits the UIR 
features in the near- and mid-IR. The near-IR spectrum is dominated by lines of H and 
He, with a hot dust continuum (Hodapp et al. 1994; Zhang & Kwok 1992). In general, 
these new data confirm the previous imaging results. In particular, the J-band image is of 
higher resolution and greater sensitivity than previous images. The halo that is detected 
in the H and A^-band images is also detected at J, and extends ~30" in diameter. One 
discrepancy between this image and the J-band image in Hora et al. is that the previous 
J image showed the east lobe significantly rotated towards the north. The location of the 
east lobe in the new J image is similar to the lobe position in the H and A^-bands, and in 
Br7 and Ha images (Louise et al. 1987). 

3.1.6. M 1-92 

M 1-92 ("Minkowski's Footprint") is a proto-typical example of a bipolar reflection 
nebula. In the visible, the fully obscured BO. 5 V (Herbig 1975) central star is seen only in 
reflected starlight scattered off two symmetrical bipolar lobes (Schmidt et al. 1978). The 
near-IR appearance is significantly different from that seen in the visible, however (Fig. 
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14b). Our near-IR images show that M 1-92 bears a striking resemblance to AFGL 618, 
which is also a bipolar proto-planetary nebula (see Fig. 13b and Latter et al. 1992). As 
in that object, M 1-92 has one bright peak situated roughly between the two lobes. The 
western lobe is brighter than the east, and the difference is much more pronounced at 
J than at K. The bright peak is not seen at visible wavelengths (Schmidt et al. 1978). 
The morphology displayed is consistent with the single-photon scattering models of Morris 
(1981) and Latter et al. (1992) as applied to AFGL 618. This implies that the polar axis of 
the nebula is highly inclined to the plane of the sky (i ~ 45°) with the western lobe being 
closest to the earth. 

3.2. Detected H 2 Emission Objects 

We imaged several objects not observed in the H 2 surveys of Zuckerman & Gatley 
(1988) and Webster et al. (1988), but no new H 2 detection was clearly made. We will, 
however, be better able to assess the H 2 content of PNe when we discuss the spectroscopic 
results in a following paper. The images of M 2-9 and AFGL 2688 acquired during the 
course of this survey have already been discussed at length elsewhere (Hora & Latter 1994; 
Latter et al. 1993) and will not be considered further here. 

3.2.1. NGC 2346 

First mapped in H 2 by Zuckerman & Gatley (1988) and recently imaged by Kastner et 
al. (1994), NGC 2346 has a remarkable butterfly-type morphology surrounding a bright 
elliptical region (Fig. 2b). The appearance of the lobes is like that of limb-brightened 
bubbles. It is likely that the H 2 emission arises from shocked gas at the interface of an 
expanding bubble of fast-moving, low-density gas interacting with the slower moving red 
giant envelope. We find the continuum to be weak, and it appears to follow the butterfly 
structure. Br7 emission is very centrally concentrated, with no detectable emission from 
the more extended elliptical region. This is different from the Ha distribution, which fills 
the central elliptical region (Balick 1987). This difference is likely, though not entirely, a 
result of insufficient sensitivity. A large amount of filamentary structure is visible in the H 2 
emission. It has been argued, based largely on the H 2 distribution, that this object and M 
57 (§3.2.2.) have very similar overall morphologies with M 57 being viewed nearly pole-on 
(Kastner et al. 1994; see also Balick et al. 1992). A contour plot of the H 2 emission is shown 
in Fig. 16. The appearance of NGC 2346 would suggest that the object is inclined by a few 
degrees such that the north lobe is pointed toward the earth. 
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3.2.2. M 57 

M 57 (NGC 6720; the Ring Nebula) has been studied in molecular hydrogen emission 
by several groups (Kastner et al. 1994; Greenhouse et al. 1988; Zuckerman & Gatley 1988). 
The main ring of this object is a strong source of H 2 emission, and has modest Br7 and 
continuum. Independent of Kastner et al. (1994), we found the emitting H 2 to extend well 
beyond the bright ring. This emission is coextensive with the visible halo, but detailed 
inspection finds that the H 2 emission is slightly offset from Ha (see Kastner et al. 1994). 
The Br7 emission primarily falls interior to the H 2 , while the continuum follows the H 2 
emission. A large amount of detailed filamentary structure is visible in our grayscale H 2 
image (Fig. 13a). These filaments point radially away from the main ring, and in at least a 
few cases seem to originate from bright knots in the ring. Such filaments might, therefore, 
reflect higher density regions created by instabilities in the flow. These data are consistent 
with the idea that M 57 is a butterfly-type nebula seen pole-on. The extent of the H 2 
emission is clearly visible in Fig. 17. Kastner et al. (1994) argue convincingly that the H 2 is 
shock excited, but we will consider this further in our forthcoming paper. A detailed study 
of the H 2 velocity structure in this object and NGC 2346 could be very illuminating. 

3.2.3. NGC 1021 

This young compact PN has a bright continuum that is largely coincident with a ring of 
Br7 emission. The near-IR morphology of NGC 7027 has been studied previously (Kastner 
et al. 1994; Graham et al. 1993b). Our H 2 image (Fig. 10a) has a strong contribution from 
He I emission, which lies in a ring inside that of the H 2 . Graham et al. (1993a) argue on 
morphological grounds that the H 2 emission must be primarily excited by UV photons, 
rather than being shock excited. Contour diagrams of continuum subtracted H 2 and Br7 
images are shown in Fig. 18. The calibration of these images agrees extremely well with 
that of Graham et al. (1993b), lending confidence to our continuum subtraction technique. 

3.2.4. J 900 

Our near-IR Br7 continuum subtracted image of J 900 (Figs. 12b and 19) is elongated 
into two lobes but otherwise featureless, much like the morphology seen in Ha and [0 ill] 
(see Schwarz et al. 1992). The continuum is elongated in the direction of the lobes. In 
contrast, we find H 2 emission to extend from the central region in almost jet-like structures. 
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At the end of the the north-west "jet" appears an arc of emission. It is likely that He I 
emission contributes to our H 2 continuum subtracted image, with a spatial distribution 
similar to that of Br7. This object has been discussed recently by Shupe et al. (1995). 

3.2.5. NGC 2440 

NGC 2440 has extended H 2 emission that does not obviously follow the structure seen 
in visible CCD images (Cutri 1994; Schwarz et al. 1992; Balick 1987). We find that the H 2 
emission extends in roughly linear structures from the center, in between the visible lobes, 
and ends in arcs that appear to almost encircle the nebula (Figs. 4a and 20). The arcs 
are not equidistant from the central region, an effect that is likely caused by orientation 
to the line of sight. Bright knots visible in our H 2 continuum subtracted image, which 
are approximately coincident with the Br7 emission, are from He 1 emission in the filter 
bandpass. 

3.2.6. AFGL 618 

AFGL 618 (Fig. 13b) is dominated by continuum emission in the visible and near-IR; 
arising primarily from scattered light emitted by the central star. It is a source of strong 
H 2 emission, but the distribution of that emission has been unknown because of the strong 
continuum (see, however, Rowlands 1993). We have successfully subtracted the continuum 
from our H 2 and Br7 images. The emitting H 2 is concentrated in the lobes (Fig. 21), similar 
to that found for AFGL 2688. No equatorial H 2 emission like that seen in AFGL 2688 is 
detected, however. How the H 2 might be excited in AFGL 618 and AFGL 2688 has been 
discussed by Hora & Latter (1994), Latter et al. (1993; 1992), and Jura & Kroto (1990). 

Br7 is confined to an unresolved, centrally located region (Fig. 21). The position and 
size of the ionized region has been discussed by Kwok & Bignell (1984). They find it to be 
~ 0'.'4 in diameter, and located roughly between the two visible lobes. This is consistent 
with our Br7 data. The integrated fluxes are /(H 2 ) = 1.1 ± 0.33 X 10~ 12 ergs s -1 cm -2 and 
/(Br7) = 1.5 ± 0.45 X 10~ 13 ergs s -1 cm -2 . These values agree very well with spectroscopic 
results (Latter et al. 1992; Kelly, Latter, k Rieke 1992; Thronson 1981). 

Broadband integrated magnitudes for these data are J = 13.7 ±0.2 mag, H = 11.4 ±0.2 
mag, and K = 8.9 ± 0.2 mag. Conditions during the broadband measurements were not 
photometric, but we consider these data to be in agreement with that of Latter et al. 
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(1992), who found a significant discrepancy with previous A^-band measurements. Our new 
J-band image is consistent with the scattering model from that study. 

3.3. Molecular Hydrogen Distribution in Planetary Nebulae 

Emission from low-lying H 2 electric quadrupole vibration-rotation transitions can be 
excited collisionally in regions with kinetic temperatures Tk ~ 1000 K. Kinetic temperatures 
like this are typically thought to be associated with the presence of moderate velocity 
shocks. For shocks to produce observable H 2 emission the velocity must fall in the range 
V s ~ 5 — 25 km s -1 , or somewhat higher for C-type shocks {V s ~ 40 — 50 km s _1 ; see 
Draine & Roberge 1982; Hollenbach & Shull 1977). Such velocities are observed in PNe, 
and H 2 excited in this way must be present at least part of the time. Ultraviolet photons 
with A > 912 A can escape from the ionized part of the nebula and then be absorbed in 
the Lyman and Werner bands of any H 2 in that region. About 10% of these absorptions 
will result in dissociation of the H 2 molecules. The remainder will lead to a cascade in 
the electronic ground state vibration-rotation levels, producing a near-IR spectrum. While 
this type of excitation can result in an easily identifiable spectrum (Black & van Dishoeck 
1987), high densities and strong UV fields can produce a spectrum that mimics that of 
shock or thermal excitation (Sternberg & Dalgarno 1989). Often without detailed spectral 
information (Hora & Latter 1994; Ramsay et al. 1993), it can be difficult to know which 
is the principal cause of the excitation. In any case, it is instructive to isolate the regions 
where H 2 is emitting in PNe. 

When present, it is evident that the emitting H 2 is always associated with the polar 
lobes of the nebulae, although it can also be present in other regions. In the case of bipolar 
reflection nebulae, such as AFGL 618 and AFGL 2688, the H 2 emission is concentrated 
within the visible lobes themselves (the equatorial emission in AFGL 2688 is discussed 
by Latter et al. 1993). In evolved butterfly-type nebulae (here we consider M 57 to be a 
member of this class), the H 2 emission comes from the outside edge of the lobes (see also, 
Kastner et al. 1994). It would appear that a transition occurs from emission within the 
lobe interior to the lobe surface, as the polar regions expand and evolve. This is the type of 
evolution that might be expected for wind-blown bubbles. 

For other objects discussed here, it is not clear how the H 2 emission is associated with 
the overall distribution of matter. While not spherically symmetric in its H 2 emission, NGC 
7027 does not obviously fall into the same class as the young bipolar reflection nebulae or 
the evolved butterfly nebulae. Its status as a very young, rapidly evolving object makes 
it likely that it should not be included in any specific evolutionary scenario of particular 
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morphologies. A similar statement can be made about BD+30°3639, which has been 
detected in H 2 elsewhere (Graham et al. 1993a; Webster et al. 1988; Zuckerman & Gatley 
1988). The cases of J 900 and NGC 2440 are even more difficult to interpret. Deeper H 2 
images are required. Visible light images of NGC 2440 show it to be a butterfly nebula 
with lobes extending in the east-west direction (Cutri 1994; Schwarz et al. 1992). The H 2 
emission is brightest between the lobes, in apparent contradiction to the trend discussed in 
the previous paragraph. In J 900 the H 2 is detected well beyond the visible manifestation 
of the nebula. The H 2 in that object extends outward in the direction to which the bright 
lobes are oriented (see Shupe et al. 1995). 

A number of studies have pointed out a strong correlation of detectable molecular 
emission with PN morphological type (e.g., this work; Kastner 1994; Huggins & Healy 1989; 
Zuckerman & Gatley 1988). From these studies, it has been made clear that butterfly or 
bipolar morphological types appear to have much larger quantities of molecular material. 
Our new data are consistent with this result - if a PN is detected in H 2 , it is a young 
bipolar reflection nebula or a more evolved butterfly nebula. As mentioned above, the 
exceptions consist of intermediate age objects that are rapidly dissociating the remnant 
AGB envelope. Note that the converse is not necessarily true - butterfly nebulae do not 
always have detectable H 2 or CO. This, however, might be more a statement of degree 
rather than contradiction. As emphasized by Huggins & Healy (1989), this correlation 
is very significant. Since the origin of nebular shapes, in particular bipolar and butterfly 
shapes, is the subject of much debate, it seems that this correlation must be investigated 
in some detail. The implications are of key importance. For example, for PNe to retain 
molecular envelopes at times significantly beyond exposure of the hot core, a minimum 
requirement is that large amounts of material must have been deposited into the molecular 
envelope. This would imply a strong correlation of mass loss rate on the AGB with 
morphological type. Also, the above mentioned molecular studies pointed out a clear trend 
for butterfly morphological types to be members of a young disk population, suggesting 
massive progenitor stars. Stars at the higher end of the AGB mass limit are capable of 
losing mass at high rates for an extended period of time before exposure of the core. It 
therefore appears that massive AGB stars produce butterfly nebulae. Thus, the close binary 
hypothesis for causing mass loss to be preferentially in the orbital plane might not be a 
necessary requirement (Zuckerman & Gatley 1988) unless all of these type of morphologies 
arise from binary systems with massive components (Han et al. 1995). 



4. Summary 



We have presented the results ol a near-IR imaging survey ol a moderately sized sample 
of planetary nebulae. Aside from the occasional presence of H 2 emission, there are no 
systematic morphological differences with PNe observed in the visible, the exceptions being 
very dusty objects like AFGL 618 and AFGL 2688. An attempt was made to make new 
detections of molecular hydrogen emission, without success. However, new images of several 
H 2 emission objects are presented. With the exception of young, compact objects, the H 2 
emission is typically found to be extended and does not always trace visible components 
of the nebulae. There is, however, some ambiguity when interpreting compact emission 
that might be confused with He I emission in the filter bandpass. Our new high-spatial 
resolution data reveal filamentary and other kinds of structures that are clearly resolved in 
many nebulae. Previous studies (e.g., Zuckerman & Gatley 1988) determined a correlation 
of H 2 emission with "butterfly" morphological type. Our results are consistent with this 
finding. It has been suggested that molecular material is only present in objects with 
massive progenitors (Huggins & Healy 1989; Zuckerman & Gatley 1988). This might well 
be the case, but the strong correlation with the morphology of PNe is still not explained, 
unless morphology is also related to progenitor mass and mass-loss rate (see Corradi & 
Schwarz 1994). The data presented in this paper will be further analyzed and combined 
with near-IR spectroscopy in a following paper. 
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Table 1 

Summary of Observations 



# 


Object 


PN G tf 1 


Type 2 


SO Filters 3 


SO Obs. dates 


UH Filters 3 


UH Obs. dates 


Fig. # 


1 


NGC 


40 


120.0+09.8 


mid. elliptical 


K, H2, Br7, cont 


15/10/92 








la 


2 


NGC 


1535 


206.4-40.5 


early round 


K, H2, Br7, cont 


15/10/92, 29/12/93 








lb 


3 


NGC 


2022 


196.6-10.9 


early elliptical 


H2, Br7, cont 


28/10/93, 29/12/93 






2a 


4 


NGC 


2346 


215.6+03.6 


mid. butterfly 


K, H2, Br7, cont 


14&15/10/92, 
10/3, 27/10 93 






2b 


5 


NGC 


2371-2 


189.1 + 19.8 


late elliptical 


H2, Br7, cont 


29/12/93 






3a 


6 


NGC 


2392 


197.8+17.3 


early round 


Ks, H2, Br7, cont 


7&9/3/93 






3b 


7 


NGC 


2440 


234.8+02.4 


late butterfly 


Ks, H2, Br7, cont 


10/3/93 






4a 


8 


NGC 


3242 


261.0+32.0 


early elliptical 


Ks, H2, Br7, cont 


10/3/93, 29/12/93 






4b 


9 


NGC 


6210 


043.1+37.7 


irregular 


Ks, H2, Br7, cont 


10/3/93 






5a 


10 


NGC 


6302 


349.5+01.0 


mid. butterfly 4 






J, K, H2, Br7, cont 


21&22/6/92, 28/6/93 


5b 


11 


NGC 


6543 


096.4+29.9 


peculiar+round halo 






H2, Br7, cont, I 


19/6/92, 28/6/93 


6a 


12 


NGC 


6572 


034.6+11.8 


late elliptical 4 ? 






K, H 2 , Br7, cont 


20/6/92, 28/6/93 


6b 


13 


NGC 


6826 


083.5+12.7 


early elliptical 
+round halo 






I 


19/6/92 


8a 


14 


NGC 


7009 


037.7-34.5 


mid. elliptical 


K, H2, cont 


15/10/92 


H2, Br7, cont 


27&28/6/93 


7 


15 


NGC 


7026 


089.0+00.3 


late elliptical 


K, H2, Br7, cont 


15/10/92, 27/10/93 


I 


19/6/92 


8 


16 


NGC 


7027 


084.9-03.4 


mid. elliptical 


K, H2, Br7, cont 


15/10/92, 27/10/93 






10a 


17 


NGC 


7662 


106.5-17.6 


early elliptical 


K, H2, Br7, cont 


14/10/92, 27/10/93 


I, Br7, cont 


28/6/93 


9 



1 



Table 1 (continued) 
Summary of Observations 



# 


Object 


PN G tf 1 


Type 2 


SO Filters 3 


SO Obs. dates 


UH Filters 3 


UH Obs. dates 


Fig. # 


18 


BD+30°3639 


064.7+05.0 


early round 






K 


20/6/92 


10b 


19 


121282+5050 




peculiar 4 






K 


20/6/92, 10/10/92 


10b 


20 


IC 289 


138.8+02.8 


mid. elliptical 


K, H2, Br7, cont 


14/10/92, 28/10/93 






11a 


21 


IC 418 


215.2-24.2 


early elliptical 


H, K, H2, Br7, cont 


15/10/92, 10/3/93 


.1 


8/2/93 


lib 


22 


IC 2149 


166.1 + 10.4 


peculiar 


H2, cont 


28/12/93 






11a 


23 


IC 4593 


025.3+40.8 


irregular 


Ks, H2, Br7, cont 


9/3/93 






12a 


24 


J 900 


194.2+02.5 


early butterfly 4 ? 


K, H2, Br7, cont 


14/10/92, 28/10, 
28&29/12/93 






12b 


25 


M 57 


063.1 + 13.9 


mid. butterfly 5 


K, H2, Br7, cont 


many 






13a 


26 


AFGL 618 


166.4-06.5 


bipolar 4 


J, H, Ks, H2, Br7, cont 


9/3/93 






13b 


27 


AFGL 2688 6 




bipolar 4 


J, H, K, H2, cont 


many 


I, H, K, H2, cont 


19&22/6/92 


14a 


28 


M 1-78 


093.5+01.4 


mid. butterfly 4 






K 


21/6/92 


15b 


29 


M 1-92 




bipolar 4 






J, K 


22/6/92 


14b 


30 


M 2-9 7 


101.8+18.0 


mid. butterfly 


K, H 2 , cont 


13/6/92, 9/3/93 


J, H, K, H 2 , Br 7 , cont 


20-22/6/92, 27/6/93 


15a 



1 Strasbourg-ESO Planetary Nebula catalog number (Acker et ah 1992). 

2 From Balick (1987) for objects which appear in both data sets (see text). 

3 Pixel scale for UH NICMOS3 data is 0'.'12 for data before 1993 June 1 and 0'.'24 after. SO data is all at O'.'6/pixel. /-band data is at O'.'08/pixel. 

4 Not in Balick (1987). 

5 Identified as a middle elliptical+round halo by Balick (1987). 

6 See Latter et al. (1993). 

7 See Hora k Latter (1994). 



2 
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Table 2 

Filter Parameters 



Instrument 


Filter 


Half-Power Wavelengths (/im) 


SO NICM0S3 Camera: 


J 


1.258 ± 0.14 




H 


1.647 ±0.17 




K 


2.205 ± 0.21 




Ks 


1.990 - 2.320 




H 2 


2.121 ± 0.01 




continuum 


2.097 ±0.01 




Br7 


2.166 ± 0.01 




continuum 


2.151 ±0.01 


UH NICM0S3 Camera: 


J 


1.238 ± 0.14 




H 


1.686 ±0.14 




K 


2.211 ± 0.20 




H 2 


2.1132 - 2.1367 




Br7 


2.1646 - 2.2053 




continuum 


2.230 - 2.290 


UH CCD Camera: 


I 


0.816 ±0.091 
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Figure 1: a) NGC 40 in A = 2.151 //m continuum, A^-band, Br7, and continuum 
subtracted Br7. H 2 is not detected. A 20" scale bar is shown. For this, and all images that 
follow, north is up and east is to the left, b) As in a) but for NGC 1535. 

Figure 2: a) NGC 2022 in H 2 , A = 2.097 //m continuum, Br7, and continuum subtracted 
Br7. H 2 is not detected. A 20" scale bar is shown, b) As in a) but for NGC 2346. The H 2 
and continuum images are displayed on a logarithmic intensity scale. H 2 is clearly detected 
in an extended butterfly-type morphology. A 40" scale bar is shown. 

Figure 3: a) NGC 2371-2 shown in H 2 , continuum subtracted H 2 , Br7, and continuum 
subtracted Br7. H 2 is not detected. A 20" scale bar is shown, b) NGC 2392 in H 2 A = 2.151 
/mi continuum, Br7, and continuum subtracted Br7. H 2 is not detected. A 20" scale bar is 
shown. 

Figure 4: a) NGC 2440 shown in H 2 , continuum subtracted H 2 , A = 2.151 //m continuum, 
and continuum subtracted Br7. H 2 is well detected. A 20" scale bar is shown for the H 2 
image, which is displayed using a logarithmic intensity scale. The 10" scale bars apply to 
the other three images. H 2 emission is found to be extended. Emission in the continuum 
subtracted H 2 image at the location of the Br7 emission is most likely He I 3 3 P° — 4 3 S 
at A = 2.114 /mi emission, which falls within the filter bandpass, b) NGC 3242 in H 2 , 
A^s-band, Br7, and continuum subtracted Br7. H 2 is not detected. A 20" scale bar is 
shown. 

Figure 5: a) NGC 6210 shown in A = 2.121 /mi continuum (H 2 is not detected), A^s-band, 
Br7, and continuum subtracted Br7. A 10" scale bar is shown, b) NGC 6302 in the UH J, 
H 2 , Br7, and continuum filters. A 20" scale bar is shown. These images are shown on a 
logarithmic intensity scale. 

Figure 6: a) NGC 6543 shown in the UH H 2 , Br7, continuum, and /-band filters. H 2 is 
not clearly detected. A 10" scale bar is shown, b) NGC 6572 in the UH K, H 2 , Br7, and 
continuum filters. A 10" scale bar is shown. 

Figure 7: a) NGC 7009 in H 2 , continuum subtracted H 2 , A = 2.151 //m continuum, and 
A^-band. H 2 is not detected. These images are shown on a logarithmic intensity scale. 
A 10" scale bar is shown, b) UH continuum and Br7 images of NGC 7009. The bottom 
frames are shown on a logarithmic intensity scale. A 10" scale bar is shown. 

Figure 8: a) NGC 6826 (lower) in the UH /-band filter and shown at different contrast 
levels. UH /-band images of NGC 7026 (upper), also at different contrast levels. These 
images are displayed on a logarithmic intensity scale. A 10" and 20" scale bar is shown for 
each object respectively, b) NGC 7026 shown in H 2 , continuum subtracted H 2 , Br7, and 
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continuum subtracted Br7. H 2 is not detected. The emission remaining in the continuum 
subtracted H 2 images is likely from He i. A 10" scale bar is shown. The Br7 images are 
shown on a logarithmic intensity scale. 

Figure 9: a) NGC 7662 shown in A = 2.121 /mi continuum (H 2 is not detected), A^-band, 
Br7, and continuum subtracted Br7. A 20" scale bar is shown. All images are displayed 
on a logarithmic intensity scale, a) NGC 7662 in UH Br7, continuum, and /-band. These 
images are shown on a logarithmic intensity scale. The /-band frames are at different 
contrast levels. A 20" scale bar is shown. 

Figure 10: a) NGC 7027 shown in H 2 , continuum subtracted H 2 , Br7, and continuum 
subtracted Br 7 . A 10" scale bar is shown, b) BD+30°3639 (upper) and IRAS 21282+5050 
in the UH A^-band. The left images are shown on a linear intensity scale, while the right 
are on a logarithmic scale. A 5" scale bar is shown for BD+30°3639 and a 5" scale bar is 
shown for IRAS 21283+5050. 

Figure 11: a) IC 289 (lower) shown in Br7 and A"-band. A 20" scale bar is shown. IC 2149 
(upper) in H 2 and continuum subtracted H 2 . The remaining emission in the continuum 
subtracted image is from He I. A 10" scale bar is shown for this object, b) IC 418 in the 
SO K, iZ, and UH J bandpasses with Br7. A 10" scale bar is shown. The continuum 
contribution to the Br7 image is moderate, and continuum subtraction (not shown) does 
not offer significantly different morphological information. 

Figure 12: a) IC 4593 shown in A = 2.121 /mi continuum (H 2 is not detected), A^-band, 
Br7, and continuum subtracted Br7. A 10" scale bar is shown, b) J 900 in H 2 , continuum 
subtracted H 2 , continuum subtracted Br7, and A = 2.151 //m continuum. H 2 is clearly 
detected. It is expected that some of the emission found in the continuum subtracted H 2 
image at the location of Br7 is from He I in the bandpass. The 20" scale bar applies to the 
top frames. The 10" scale bar applies to the bottom frames. 

Figure 13: a) M 57 in H 2 , A = 2.151 //m continuum, Br7, and continuum subtracted Br7. 
H 2 is well detected. A 40" scale bar is shown. The H 2 image is on a logarithmic intensity 
scale, b) AFGL 618 shown in continuum subtracted H 2 , continuum subtracted Br7, and in 
the J and K bandpasses. A 5" scale bar is shown. The bottom two frames are shown on a 
logarithmic intensity scale. 

Figure 14: a) AFGL 2688 in H 2 , continuum subtracted H 2 , A = 2.151 //m continuum, and 
in the UH A^-band. The reduced 10" scale bar applies to the A^-band frame only. These 
data are discussed in detail elsewhere (Latter et al. 1993). b) M 1-92 shown in the UH J- 
and A^-bands. The upper images are shown on a logarithmic intensity scale. The lower are 
on a linear scale. A 5" scale bar is shown. 
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Figure 15: a) M 2-9 in the UH H 2 , Br7, continuum, and A^-band niters. A 10" scale bar 
is shown. A detailed discussion can be found in Hora & Latter (1994). b) M 1-78 in the 
UH A^-band. the intensity scale is logarithmic, with the frames shown at different contrast 
levels. A 10" scale bar is shown. 

Figure 16: Contour plot of NGC 2346 in the SO H 2 filter. Contours are from 5.6 X 10 -16 
to 8.3 X 10 -15 erg cm -2 s -1 arcsec -2 with an interval of 5.6 X 10 -16 erg cm -2 s -1 arcsec -2 . 
The continuum is weak, and has not been subtracted from this image. 

Figure 17: Contour plot of M 57 in the SO H 2 filter. Contours are from 1.9 X 10 -16 to 

5.3 X 10 -15 erg cm -2 s -1 arcsec -2 with an interval of 5.6 X 10 -16 erg cm -2 s -1 arcsec -2 . 
The continuum is weak, and has not been subtracted from this image. 

Figure 18: a) Contour plot of NGC 7027 in the SO H 2 filter. Contours are from 2.5 X 10 -15 
to 1.4 X 10 -13 erg cm -2 s -1 arcsec -2 with an interval of 4.2 X 10 -15 erg cm -2 s -1 arcsec -2 . 
The bright inner ring is from He I 3 3 P° — 4 3 S at A = 2.114 /mi emission, which falls 
within the filter bandpass. However, the more extended emission is H 2 . b) As in a) but 
for the SO Br7 filter. Contours are from 1.1 X 10 -14 to 1.0 X 10 -12 erg cm -2 s -1 arcsec -2 
with an interval of 2.8 X 10 -14 erg cm -2 s -1 arcsec -2 . These images have been continuum 
subtracted. 

Figure 19: Contour plot of J 900 in the SO H 2 filter. Contours are from 1.4 X 10 -16 to 

1.4 X 10 -15 erg cm -2 s -1 arcsec -2 with an interval of 1.1 X 10 -16 erg cm -2 s -1 arcsec -2 . 
The continuum is compact and has not been subtracted from this image. 

Figure 20: Contour plot of NGC 2440 in the SO H 2 filter. Contours are from 4.2 X 10 -16 
to 2.6 X 10 -15 erg cm -2 s -1 arcsec -2 with an interval of 2.8 X 10 -16 erg cm -2 s -1 arcsec -2 . 
The continuum is not extended with the H 2 emission, and has not been subtracted from 
this image. 

Figure 21: a) Contour plot of AFGL 618 in the SO H 2 filter with continuum subtracted. 
Contours are from 2.8 X 10 -15 to 2.2 X 10 -13 erg cm -2 s -1 arcsec -2 with an interval of 
8.3 X 10 -15 erg cm -2 s -1 arcsec -2 . b) As in a) for the SO Br7 filter. Contours are from 
1.9 X 10 -15 to 2.2 X 10 -13 erg cm -2 s -1 arcsec -2 with an interval of 8.3 X 10 -15 erg cm -2 
s -1 arcsec -2 . 



